Diffusion in infinite and semi-infinite lattices with long-range coupling 
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We prove that for a one-dimensional infinite lattice, with long-range coupling among sites, the 
diffusion of an initial delta-like pulse in the bulk, is ballistic at all times. We obtain a closed- 
form expression for the mean square displacement (MSD) as a function of time, and show some 
cases including finite range coupling, exponentially decreasing coupling and power-law decreasing 
coupling. For the case of an initial excitation at the edge of the lattice, we find an approximate 
expression for the MSD that predicts ballistic behavior at long times, in agreement with numerical 
results. 
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The physics of discrete systems have been a topic of 
interest for many years, because it can give rise to com- 
pletely different phenomenology in comparison with that 
present in homogeneous continuous systems. In partic- 
ular, discrete periodic systems are found in many differ- 
ent contexts including condensed matter physics, optics, 
Bose-Einstein condensates and magnetic metamaterials 
among others. Under the appropriate approximation, 
they can all be described by some variant of the dis- 
crete Schrodinger (DS) equation pQ. In that way, many 
of these systems display the phenomenology common to 
periodic systems such as the presence of a band structure, 
discrete diffraction, Bloch oscillations, dynamic localiza- 
tion, Zcncr tunneling, to name a few. 

Usually, the DS equation is used in the weak-coupling 
limit, where the interaction among sites includes nearest- 
neighbors only. This is a good approximation in cases 
where the coupling among sites decays very quickly 
with distance, like the exponentially-decreasing coupling 
found in optical waveguide arrays. However, there are 
cases where it is advisable to go beyond this approxima- 
tion. An example of that is a split-ring resonator (SRR) 
array, where the interaction among the basic units is 
dipolar in nature and therefore, the coupling decreases 
as the inverse cubic power of the mutual distance. 

When coupling beyond nearest-neighbors are consid- 
ered, the number of possible routes of energy exchange 
increases. The dynamical evolution of excited pulses in fi- 
nite ID and 2D lattices with anisotropic couplings and up 
to second nearest-neighbor couplings, has been explored 
in [5] by means of the Green function formalism. Ex- 
perimental observation of the influence of second order 
coupling in linear and nonlinear optical zig-zag waveg- 
uide arrays has been recently carried out [3J. In a dif- 
ferent context, a recent work 0], shows that long-range 
coupling in low dimensional system can induce a phase 
transition from delocalized to localized modes. Another 
interesting scenario that can be modeled as a discrete 
system with long-range couplings is that of complex net- 
works [5 , where the distances between nodes are not nec- 
essarily physical. 

In this Letter we carry out an analytical and numerical 
study on the diffusion of an initially localized pulse prop- 



agating in a one-dimensional discrete periodic lattice, in 
the presence of arbitrary long-range couplings. We focus 
on two cases of interest: A delta-like pulse in the bulk, 
and a delta-like pulse at the edge of the lattice. The ana- 
lytical work centers on the evaluation of the mean square 
displacement of the excitation, which is obtained in exact 
form for the bulk excitation, and in an approximate form 
for the edge excitation. 

Let us consider the (dimensionless) discrete 
Schrodinger(DS) equation, describing the evolution 
of an excitation along a one-dimensional periodic lattice: 

% ~([ z + v n,m u m = (1) 

where u n is the complex amplitude of the excitation at 
the nth site, z is the evolution coordinate ('time' in the 
tight-binding model for electrons, or 'longitudinal dis- 
tance' for coupled waveguide arrays in optics ). Matrix 
element V n , m denotes the coupling between the nth and 
mth sites. This matrix is periodic in space and obeys 
V n ,m = V m ,n = Vj„_ m |. Eq. |l]) conserves the norm 
P = X^oo l u »( z )| 2 which can then be set, without loss 
of generality, as unity, P = 1. The dispersion relation of 
the linear waves is obtained by inserting a solution of the 
form u n = Ae l ( kn+Xz ^ in Eq. (1), obtaining 

oo 

A = Vn, m e lk{rn - n) = 2 J2 V m cos(mk). (2) 

The convergence of this series for all k values constrains 
V n to decrease faster than 1/n. From Eq.pl) we im- 
mediately obtain some basic properties of X(k): A(fc) = 
A(-fc), A(fc) = A(fc + 27r<7), q G Z, and <9 fc A(fc)|fc=o = 
d k X(k)\ k=±7T = 0. 

The dynamical evolution of an initially localized pulse 
in a lattice, u n (0) — AoS njTlo , can be monitored through 
the mean square displacement (MSD) of the excitation: 

oo oo 

<n 2 )^n 2 K(z)| 2 /£M z )! 2 ' (3) 

M M 

where |no| S> and M = — oo for an infinite lattice, or 
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FIG. 1: Discrete diffraction pattern for a bulk spread of a 
delta-like initial condition at Vz — 20 for different values of 
dispersion parameter a in a lattice with coupling exponen- 
tially decreasing. 



no = and M = for a semi-infinite lattice. 

Diffusion in the bulk. In this case u n (0) = AoS n _ no , 
where n is far away from the boundaries of the lat- 
tice. By combining Fouricr-and-Laplace transforms to 
Eq. (1), followed by the corresponding transformation 
back to space n and z coordinates, one obtains a formal 
expression for u n (z) 



u n (z) 



2tt 



i(kn-X(k)z) ^ 



(4) 



where X(k) is the dispersion relation. Next we proceed 
to insert Eq.Q into Eq.([3]), and after using the general 
properties of A& , one obtains after some algebra: 



(n 2 ) 



2W-A dk J 



(5) 



Equation |5]) implies that the propagation of the excita- 
tion is ballistic at all times, with a 'speed' that depends 
upon the smoothness of the dispersion relation. This re- 
sult is also valid in any dimension d, where it can be 
easily proven that 



<n 2 ) = 



1 



(V fe A(fc))Vfc 



FBZ 



(6) 



where the integral is taken over the first Brillouin zone 
(FBZ), with volume v. An equivalent expression to Eq. 
^ in "real space" is obtained by replacing ^ in ^ 



<n 2 > 



(7) 



where m is the relative position if a lattice node from an 
arbitrary site n taken here as without loss of generality 
because of periodicity, and Vo,m is the coupling between 
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FIG. 2: Speed of ballistic propagation of initially localized 
excitation, as a function of dispersion parameter. The solid 
(dashed) curve corresponds to exponential (power-law) cou- 
pling case. For power-law coupling g(a) = 2£(2(a — 1)). 



sites and m, where V n m is periodic in space and obeys: 

Let us compute now in detail the MSD for several cases 
of interest. 

(a) Second-order coupling. In this case, V\ = V, V% = 
j3V and K;>2 = 0, and the dispersion relation is A = 
2V(cos(fc) +/3cos(2£;)). This implies, according to Eq.|5l) 



(n 2 ) =2{l + Ap 2 ){Vz) 2 



(8) 



Thus, the inclusion of a second-order coupling (only) al- 
ways increases the speed of the transversal diffusion. 

(b) Exponentially- decreasing coupling. In this case, 
V n . m = I/e~ a (' Tl ~ m ' _1 ' ) , where V is the coupling between 
nearest-neighbors sites and a is the long-range parame- 
ter. The linear dispersion relation is 



A = V 



e a cos(fc) — 1 
cosh(a) — cos(fc) 



(9) 



Figure [T] shows some snapshots of the spatial transver- 
sal profiles at a given longitudinal propagation distance 
z, for different values of the dispersion parameter. We 
note that, at small a values, a fraction of the initial ex- 
citation seems to 'linger' at the initial positioned, m 
a short of quasi-selftrapping, while the "untrapped" por- 
tion propagates away from the initial site faster than in 
the nearest-neighbor case. These observations can be put 
on a more rigorous basis by computing the mean square 
displacement directly from Eq. |5|, using Eq. (|9|. We 
obtain: (n 2 ^ = g(a)(Vz) 2 where, 



g(a) = — coth(a)(coth(a) + l) z 



(10) 



This implies a ballistic speed g(a) > 2, for all finite a. 
However, it only becomes noticeable larger than two for 
a < 2 (Fig. |J). 

(c) Power-law coupling. Another popular case of long- 
range interaction is the power-law coupling V n}Tn = 
V/\n— m\ a . In this case the dispersion relation is given in 



3 



terms of a polylogarithm function X(k, a) = V (Li, 

+Li a (e~ ik )) = 2VY,2=i cos i mk )/ ma - For a = 
can explicitly written as 

A(M) = -^log|2-2cos(fc)| 



,(e ife ) 
1 we 



(11) 



with the logarithm diverges for k = 0. On the other 
hand, for a > 1 the polylogarithm remains bounded at 
|fc| < 7T. Moreover, in this case we can use Eq. to 
easily calculate the mean square displacement, obtaining 



(n 2 )=2C(2(a-l))(^z) 2 



(12) 



where £ is the Riemann zeta function. Figure [2] shows the 
"speed" of diffusion as a function of a and compares it 
with the "speed" obtained for lattices with exponentially- 
decreasing coupling. The discrete diffraction pattern 
(Fig.l) observed in case (b), is also observed in this case. 

In all three cases examined, the speed of diffusion 
is greater than in the case with coupling to nearest- 
neighbors only. For the exponential and power-law cases, 
we note that as the range of the interaction is increased, 
the pulse experiences a sort of quasi-localization at the 
initial site. At the same time the speed g(a) increases 
above 2 and tends to diverge at a = (Fig. 2). We can 
understand this quasi-localization phenomenon by ana- 
lyzing the dispersion relations more closely. In cases (b) 
and (c), as the range increases, the dispersion gets flatter 
and flatter in the vicinity of = ±tt, signaling the emer- 
gence of linear modes with very small group velocities. 
As a is increased further, this vicinity grows and most of 
the linear modes acquire negligible velocity, with the ex- 
ception of a small vicinity of k = 0, where the concavity is 
very high, giving rise to long wavelength modes with high 
velocity. Since our initial delta-like condition is a super- 
position of all these linear waves, the pseudo-localization 
phenomena can be understood as due to those modes 
with zero velocity, while those few and fast long wave- 
length modes, give rise to the wings that escape at high 
speed. 

In the formal limit a — and for a finite number of 
sites N, all sites are connected to one another, and the 
high degree of degeneracy forces the localization of the 
wave function through the execution of incomplete oscil- 
lations between the initial excited site and all the other 
ones. In this case the effective dynamics can be mapped 
to the dynamics of an asymmetric dimer[B]. This system 
has been solved in closed form and shows partial linear 
localization at the initial site for finite N, which becomes 
complete in the limit N — > oo [7] . 

Diffusion at the boundary. In this case u„(0) = Ao5 n ^, 
i.e., the initial excitation is at the very edge of the lat- 
tice. This system can be viewed as an infinite lattice with 
boundary conditions: u n (z) = for n < 0, for all z. The 
lack of translational symmetry prevent us from writing 
a general expression for u(z), similar to Eq.([5]). How- 
ever, the case of nearest-neighbors coupling only, can be 
solved in closed form using the method of images [5]. The 
boundary conditions are obeyed if we placed an opposite 
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FIG. 3: Discrete diffraction pattern for a surface spread of 
a delta-like initial condition at Vz = 20 for different val- 
ues of dispersion parameter a in a lattice with exponentially- 
decreasing coupling. The solid line corresponds to approxi- 
mation ( 13 1. 



delta-like source at n = —2. Then, the propagation for 
ri > is then the superposition of the evolution of both 
sources 



u n (z) 
An 
2tt 



( (l-e- 2lk )e l{kn - x{ - k)z) dk, (13) 

J — IT 



where index s denotes a semi-infinite array. However, 
When long-range couplings are considered, it is no longer 
possible to satisfy the boundary conditions by the images 
method. The reason has to do with the fact that, even 
though we have U-i(z) = for all z, the fields of each 
semi-infinite arrays are still coupled by the range of the 
interaction. 



However, by comparing the predictions of Eq. ( 13 1 with 



numerical results for the dynamical evolution, we find 



that Eq.(13) still provides a good approximation to the 



general case, by just substituting the corresponding dis- 
persion relation for a given long-range coupling. After 
using general properties of the dispersion relation, one 
obtains an asymptotic expression 



2\ s 



<„') 



(14) 



valid for z ^> 1. 

As expected, the asymptotic behavior is ballistic in all 
cases, after some transient time, related to the repulsive 
effect of the boundary. Equation pT] ) can also be ex- 
pressed in real space, by replacing Eq.([7]) in Eq 
obtaining 



(141, 



2\ s 



+ m 2 V^ - 2 m(m + 2)V m V m+2 

(15) 
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FIG. 4: MSD in logarithmic scale for exponentially-decreasing 
coupling. Black line shows the approximate solution obtained 
using the method of images, Eq.(13l, while the dashed line 
shows the slope associated to z 2 (to guide the eyes only). 
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FIG. 5: Square of speed of ballistic propagation of initially lo- 
calized excitation, as a function of dispersion parameter a for 
the exponentially-decreasing coupling. Triangles correspond 
to numerical values, continuous line to approximation (161 
and dashed line to the exact bulk solution, Eq. <J 1 1 



This results are confirmed by direct computations of 
cases (a), (b) and (c), using 'exact' Eq.(msdsur) with 
an edge excitation as an initial condition, and compar- 
ing it with the asymptotic formulas. Fig|4] shows this 
comparison of the MSD for the exponentially-decreasing 
coupling case. In this case, the asymptotic behavior is 



(n 2 ) s 3coth(a)(Vz) 2 



(16) 



In Fig. (J5j) we show a comparison between this approxi- 
mation and the numerical results. We see that for a > 1 
there is a good correspondence between the numerical 
and analytical approach, while for a < 1, the exact so- 
lution for the bulk, Eq. (10 1, constitutes a better fit 
than the surface approximation, Eq.(16). This is under- 
standable because in the long-range coupling regimen, 
the boundary loses its meaning as such and the edge site 
behaves effectively like a bulk site. 



In conclusion, we have examined the propagation of an 
excitation initially placed in the bulk and at the bound- 
ary of a one-dimensional lattice with long-range coupling. 
For the bulk case, we find a closed-form expression show- 
ing that the propagation is ballistics at all times, while 
for the edge excitation, we obtain an asymptotic approx- 
imation based on the method of images, that predicts 
ballistic propagation at long times, in agreement with 
numerical results. 
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